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The physical properties of a supersaturated binary solution such as its density p, shear viscosity 77, and 
solute mass diffusivity D are dependent on the solute concentration c: p=p(c), 7/ = r/(c), and D =D{c). 

The diffusion boundary layer equations related to crystal growth from solution are derived for the case 
of natural convection with a solution density, a shear viscosity, and a solute diffusivity that are all depen- 
dent on solute concentration. The solution of these equations has demonstrated the following, (a) At the 
vicinity of the saturation concentration c 5 the solution shear viscosity 77 depends on p as 
Vs=V(ps ) cc p W2 (c s ). This theoretically derived result has been verified in experiments with several aque- 
ous solutions of inorganic and organic salts, (b) The maximum solute mass transfer towards the growing 
crystal surface can be achieved for values of c where the ratio of d \n[D{c)/dc] to d \n[rj(c)/dc ] is a 
maximum. 

PACS number(s): 61.30.-v, 64.60.Cn, 64.60.My, 47.15.Cb 


I. INTRODUCTION 

The problem of crystal growth from supersaturated 
solutions is well known. In recent years there has been 
considerable interest in crystal growth under microgravi- 
ty conditions. However, there is no appropriate formal- 
ism that accounts for both the hydrodynamic and the 
thermodynamic aspects of the problem. The develop- 
ment of such a formalism is of paramount importance in 
understanding how the static and the dynamic charac- 
teristics of a solution, such as its density, viscosity, and 
diffusivity, are related to each other at concentrations at 
saturation and above. 

Current analytical approaches for the descriptions of 
crystal growth from supersaturated solutions (see [1-3] 
and references therein) neglect most of the significant 
features of the solution’s metastable state. These include 
the nontrivial dependence of the solution density p, the 
shear 77 and the bulk f viscosities, and the solute mass 
diffusivity D on the solute concentration c. The usual 
practice [1-3] has been to assume that these physical 
characteristics are constants that are independent of 
solute concentration. However, recent studies [4-11] of 
supersaturated aqueous solutions of inorganic and organ- 
ic salts have demonstrated that the solution density p, the 
shear viscosity 17, and the solute mass diffusivity D have a 
nontrivial dependence on the solute concentration c: 
p=p(c), rj = j](c) y and Z> =Z) (c), which becomes more 
significant with deeper penetration into the metastable re- 
gion. For example, the diffusivity D(c) declines to zero 
at the spinodal line that separates metastable and unsta- 
ble states. These facts have necessitated the development 
of an appropriate formalism to describe crystal growth 
from supersaturated solutions taking into account the 
dependence of these physical properties of supersaturated 
solutions on solute concentration. 
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II. NATURAL CONVECTION EQUATIONS 
AND SOLUTION METASTABILITY 

Let us consider the situation where a supersaturated 
solution is mixed as a result of natural convection, which 
arises due to a depletion of solute concentration near the 
growing crystal surface. This depletion results in a 
change in solution density and leads to the appearance of 
concentration flows. However, because the solution has a 
finite viscosity and sticks to the crystal solid interface, an 
unmixed boundary layer is formed. Within this layer, ad- 
jacent to the growing crystal surface, the solution can be 
assumed to be stationary and the solute mass transfer is 
achieved only by means of ordinary diffusion. In this pa- 
per we assume that the crystal surface and the solution 
have the same temperature and thus there is no heat 
transfer. 

Natural convection arises only when there is a change 
of solution density occurring in a gravitational field and 
only in those cases where the density gradient formed is 
perpendicular to a gravitational field or when the solution 
density increases from the bottom upward. The magni- 
tude and the distribution of hydrodynamic and 
diffusional flows depend primarily on geometry and, in 
particular, on the shape and orientation of the growing 
crystal surface. Generally, on the basis of the Bjerkness 
theorem [12], one may expect that natural convection 
will occur in such a way that surfaces of equal solution 
density are oriented perpendicular to surfaces of equal 
pressure. In this paper we consider the simplest case 
where the growing crystal surface is a smooth vertical 
plate placed in a gravitational field. The case of horizon- 
tal orientation of the growing surface is considerably 
more complicated and will be studied later. 

Solution metastability implies that its density, shear 
and bulk viscosities, and solute mass diffusivity are func- 
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tions of solute concentration. Due to the crystal growth 
process the solute concentration is different far from the 
growing surface than at the crystal-liquid interface. It is 
usually assumed that an entire change of solute concen- 
tration occurs within the diffusion boundary layer (DBL) 
8 dlff . At the growing crystal surface the solute concentra- 
tion c must be greater than or equal to the saturation 
concentration c s (c>c s ) at the system temperature and 
pressure. Outside the boundary layer the solute concen- 
tration is a constant (c=c 00 ), although it fluctuates in 
time and space. Thus, within the boundary layer 8 diff , for 
the case when the solution density pic), the shear 77 (c) 
and the bulk f(c) viscosities, and the solute diffusivity 
Die) are weak functions of the solute concentration c one 
can write the following Taylor expansions: 


p{c)=p{c x ) + 


3p(c) 

3c 



Ac[l + 0(1)] , 


( 1 ) 


7](c) = 7]{C x ) + 


377 (c) 

3c 


Ac[l + 0(1)] , 

c — c „ 


( 2 ) 


D(c) = D(c x ) + 


3£>(c) 

3c 


Ac[ 1 + 0( 1 )] , 

C = C 


(3) 


where Ac=c — c x . The weakness concept for any 
infinitely times differentiable function /(c) assumes that 


— In— L^£L A c «n where n > 2 . 

dc dc n ~ ] c=c x 


In spite of the possible types of dependence of the func- 
tion /(c) on c, fulfillment of the weakness concept is al- 
ways expected for the case of small supersaturations 

i I max I C 00 i • 

Let us take x 2 =0 as the surface (010) of the growing 
crystal plate, directing the x 2 axis into a solution and the 
x x axis vertically upward. The lower ledge of the plate 
corresponds to the value x,=0. As a supersaturated 
solution flows parallel to the crystal plate, the thickness 
of the boundary layer 8 diff increases with the distance x { 
from the lower ledge where flow meets the plate. There- 
fore, the thickness 8 dlff of the DBL can be considered as 
an increasing function of x x : S d jff— 5 dlff (x , ). Thus the 
general Navier-Stokes equations describing the gravity- 
induced two-dimensional steady laminar motion of 
an incompressible supersaturated solution within 
[x 2 < Suffix x )] and beyond [x 2 > 8 di(! <x 1 )] the DBL ac- 
quire the form 


2 3 n,.- 

2~=s,,g[p(c)-p(0]m diff u 1 >-* 2 ) , 

j — 1 0X J 

3c, 3c 2 

— L + — -=0 , 

3*! 3x 2 


(4a) 

(4b) 


where 


solution momentum flow, p is the solution pressure, Uj 
and v 2 are the solution velocity components, g is the 
gravity acceleration, and 6{6 d ^x j ) — x 2 ) is the unit step 
function (Heaviside function) equal to 1 for )>x 2 

and equal to 0 for 8 diff (x , )<x 2 . Equation (4a) assumes 
that solution is isotropic, i.e., 77 (c) is the scalar quantity. 
The continuity equation (4b) assumes that a supersaturat- 
ed solution is an incompressible fluid. Utilizing the fact 
that the DBL thickness S diff ( a: 1 ) is very small compared 
to the length L along a crystal plate, one can significantly 
simplify Eq. (4a). Thus, in the particular case where 
8 din<*i)/ L « 1 , the equation describing solution motion 
along the x { axis acquires the form 


pic) 


3c, 3c, 

+ v 2 


3x 


dx 1 



3c, 


+g[p(c)-p(c oo )]0(8 di ^x I )-x 2 ) ■ (5) 


This form implies the following assumptions imposed on 
solution motion: (a) viscosity exerts a significant 

influence only within a boundary layer, (b) flow velocity 
must become zero at the solution-crystal interface, (c) the 
solution retardation in a boundary layer is caused by 
viscous forces alone, and (d) the pressure change in a 
boundary layer is determined by its change outside and in 
the natural convection case is negligible since solution 
outside a boundary layer is almost stagnant. Substituting 
into Eq. (5) expansions (1) and (2), which account for the 
dependences of the solution density and the shear viscosi- 
ty on solute concentration and keeping only the main 
terms proportional to Ac and 3( Ac )/3x 2 one can rewrite 
Eq. (5) as 


3c 1 
3x, 


-Vt + 


3c j 
3 2 Cj 


3x] 




3 2 c t 

3x? 


Ac#(8 diff (x j )— x 2 ) 


" V ^3^ 


3cj 


3xn 


Ac6>(8 difr (jc , )-x 2 


( 6 ) 


where 


’ 


_ 31n[p(c)] 
3c 


n ij = —pb ij +p(c)v i Vj+T](c) 


3K, dvj 

3x j 3x, 


*7 = 1,2 . 


_ 31n[77<c)] 


3c 


In the above expressions II i7 is the density tensor of the Boundary conditions for Eq. (6) have the form 



52 


CONCENTRATION DEPENDENCE OF SOLUTION SHEAR . . . 


807 


III. SOLUTION OF THE DERIVED EQUATIONS 
FOR THE CASE OF NATURAL CONVECTION 
IN SUPERSATURATED SOLUTIONS 



(7b) 


Condition (7a) reflects the fact that solution sticks to the 
crystal solid interface, whereas condition (7b) describes 
solution stagnancy at the infinite distance from the DBL 
edge. Thus Eqs. (4b) and (6), together with boundary 
conditions (7a) and (7b) describe the hydrodynamic as- 
pects of natural convection corrected by accounting for 
solution metastability in the system supersaturated solu- 
tion plus growing crystal surface . 


In order to simultaneously solve Eqs. (4b) and (6) let us 
introduce the stream function Wix x ,z) = G { x\ /A f(z) and 
z — (S 2 x f 1/4 x 2 defined by the relations [12,13] 




»2 


d y V{x ]y z ) 
0*2 


=e x e 2 x\ n f(z) 


0*P(x j ,Z ) 

dx { 


^-V /4 [3 fU)-f'(z)) , 


( 8 ) 

(9) 


where the dimension constants Q x and (S 2 and the func- 
tion f(z) are the quantities to be defined later. The intro- 
duction of the stream function 'Pfxpz) allows us to 
rewrite Eq. (6) in the form 


j[/'<* >] 2 -4/< z )/"(*)= -v„ ^/'"<Z)[ 1 ■ +u. )0<8 diff U , )~x 2 )Ac ] 


<?, 


e 




2 r „. , 9Ac , oo . 
Voodoo e f (Z) 3z + ((?](?2) 2 AC 


^(S di ff<jc , )— JC 2 ) • 


( 10 ) 


The equation obtained above can be satisfied only when the solute concentration c is understood as a function of the 
variable z: c=c(z). In addition, as described above, there exists the equality c(z)=c co for the case when x 2 > b^Xj ). 
This means that along the DBL edge the variable z is constant: z ~z x . Therefore, after accomplishing straightforward 
but cumbersome calculations we come to the following important conclusions: (a) the DBL thickness S difr (x ] ) is defined 
by the relation 




1/4 


(11a) 


and (b) in the vicinity of the DBL edge there exists the following expansion of A c =Ac(z): 

dciz) 


Ac(z) = rjz-z o0 )9(z x -znO((z-z x ) 1 ) , r oc = 


3z 


(lib) 


where diz^—z) is the Heaviside function [see definitions for expression (4a)]. Substituting now expansion (lib) into 
Eq. (6), we obtain the nonlinear differential equation for the function /(z), 


\[f , (z)] 2 -~f(z)f , '(z)=-v x -^f , ' t (z)[\+(co Vi -K oc )rjz-z a0 )fl(z M -z)] 

<?2 g^ocXoc 

+ vc # ~r al /"(z)9(z fi0 —z)— (z-z^ )0(z x> -z) . 




(e x e 2 ) 2 


( 12 ) 


Let us look for a solution of this equation in the form 

f(z)=fJz)G(z O0 -z)+f b (z)e( Z -z QO ) y (13) 

where the functions f w (z) and f b (z) give solutions of Eq. 
(12) within and beyond the DBL. Equations for the func- 
tion f w (z) and f b (z) have the forms 


|[/ft(z)] 2 -iA(z)//(z)=-v 00 |^-/;'(z) . (14b) 

Boundary conditions for these equations given in terms of 
the functions f w (z) and f b (z) can be derived from rela- 
tions (7)-(9) in the form 




/u,< z )L=o = /»( z )l*=o=0 . 


(15a) 


-v x ^r;’(z)[\+(<o x -Kjrjz-z x )} 


+v 00 K 00 |^-r a0 /;:(z)- 


«? i e 2 ) 2 


(z-z.) , (14a) 


f b (z)\ z = x =f' b (z)\ z=x =0 . (15b) 

One of the possible solutions of Eq. (14a) subjected to 
the boundary condition (15a), can be found as the series 
expansion 
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/.<*>= i . 

«=o n * 

where 

/o=/.=/ 2=0 , 


/j — 


v.<?,<?|[i+(ic.-fl).)r ce z. 


U-~rr. ^ r /3 > 


/s — 


n-u^-cOr^i 
z^-z^ 

1 +(*„-<», ir**. 


A , 


A 


Ak^ — 2xd n 


i +u w -co x ir^z,, 


/s > 


(16) 


_ (n — 1 )*<* — (a — 2)^ 
fn +l ~ l+u.-wjr.z, /n 

- V Ju ~ 2)! v, [5w (» -m)-l«(H - 1 )] 

Xfmfn-m > 


where « > 6. This solution can be considerably simplified 
if one assumes that the solution constant z ^ can be 
defined as = l/(r oo #c 00 ). Such an assumption leads to 
the vanishing of all expansion coefficients f n (n >4) leav- 
ing only one nonzero term in series (16) corresponding to 
n =3, 


/uA)=i/> 3 , fi = 


go* 

(2k x -(0 0 


(17a) 


e,=- 


g<»o 


72v 2 A(o x -2k x ) 


1/4 


(18) 


It is noteworthy that further agreement (patching) of the 
solutions f w (z) and f b (z) is impossible since their first 
derivative and all consequent derivatives are not equal 
along the DBL edge. Therefore, one may conclude that 
the problem general solution /(z), given by relation (13), 
has discontinuity points along the DBL edge. To explain 
this conclusion let us find velocity components v { and v 2 
within and beyond the DBL. For this purpose results 
(17a) and (17b) for the functions f w (z) and f b (z) have to 
be substituted into the corresponding expressions (8) and 
(9) for the velocity components: 




g^QO 2 

2 vJ2k x -o) k ) X2 


v 2 =0 (within) , 


x \ 1 

V\- 12v w — , u 2 = 12v oc — (beyond). 
x\ x 2 

It follows from the expressions obtained for the velocity 
components v { and v 2 that within the DBL the tangent 
velocity component v ] is infinitely times greater than its 
normal component v 2 , whereas beyond the DBL these 
components become of the same order. This corresponds 
to the following well-understood phenomenon of the 
DBL separation: (a) within the DBL the solution motion 
occurs along the crystal surface without side deviations 
and (b) along the DBL edge the tangent and normal ve- 
locity components become of the same order and the 
DBL separation takes place. Thus the DBL edge can be 
identified with the discontinuity line in solution motion 
[13,14]. 


One of the possible solutions of Eq. (14b) subjected to 
the boundary condition (15b) is 

^2 1 

f b (z)=- 12v . (17b) 

L j Z 

Both solutions f w (z) and f b (z) obtained correspond to 
connective flows within and beyond the DBL and should 
coincide along this layer edge )=f b (z ao ). This 

condition provides a way to determine the constant (? 2 in 
the form 


IV. THE DBL THICKNESS AND RELATIONSHIP 
BETWEEN SOLUTION SHEAR VISCOSITY 
AND DENSITY AT SATURATION POINT 


Substituting relation (8) into expression (11a) for the 
DBL thickness ) we obtain the result 


W*i) = 


72v 2 J2k x -co 00 )x 1 


l*l< 


1/4 


(19) 


In the expression above we have taken into account that 


TABLE I. Dependence of solution density (g/cm 3 ) on solute concentration (mass fraction). The sat- 
uration point is taken at 25 °C and normal pressure. 


Aqueous 

solution 

(? of 

^7 sample^ <* ) 

a 0 


Saturation point 
at 25 °C 

and normal pressure 
(in mass fraction) 

Intersection point 
between 17 ( 0 *) 
and r/ samplc (c 00 ) 
(in mass fraction) 

KC1 

0.904 

1.013 

0.571 

0.264 

0.268 

ADP 

1.739 

0.998 

0.564 

0.283 

0.253 

KDP 

1.294 

0.996 

0.744 

0.200 

0.195 

TGS 

1.437 

1.005 

0.367 

0.231 

0.215 

glycine 

1.329 

1.000 

0.343 

0.198 

0.183 
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the x | axis and the gravity acceleration g are oppositely 
directed. Thus, by means of the result (19), we have de- 
rived how the hydrodynamics of natural convection are 
related to thermodynamic metastability of supersaturated 
solutions. In expression (19) metastability effects are tak- 
en into account through the dependence of the solution 
density pic ^ ) and the shear viscosity r\(c ^ ) on the bulk 
solute concentration c x . 

An analysis of expression (19) obtained for the DBL 
thickness S^Xj) allows the following conclusions. 
First, we have obtained the well-known result that 
£>diff^*i) grows as ix ] /\g\) W4 (see [12-14] and references 
therein). Second, we have derived how the DBL thick- 
ness S^Xj) depends on the bulk solute concentration 
c ^ via such solution static and dynamic characteristics as 
its density pic O0 ) and viscosity J]ic ^ ). This dependence 
provides an opportunity to relate to each other the super- 
saturated solution static and dynamic characteristics. In 
particular, expression (19) allows a relationship between 
the density pic s ) and the viscosity 7]{c s ) of supersaturated 


1.8 

1.7 

1.6 

1.5 


; y]( c yj Mnpit {c) ( cP ) 


(a) 


1.2 
1.1 \ 

0.08 


TGS 


saturation point 


20 % 


solute concentration , c (in mass fraction) 


r](c),7i <c)( c p) 

' sample 


(b) GLYCINE 


solutions. It is apparent that a boundary layer should 
vanish at the saturation point since at this point solution 
and the growing crystal surface are in thermodynamic 
equilibrium. As it follows from the result (19) obtained 
for bji^Xj ), such a situation is possible only when the fol- 
lowing equality is satisfied: 

=c = 2*Jc =c • ( 20 ) 

oo s CD s 

An analysis of Eq. (20) gives that at saturation point the 
solution density pic s ) and viscosity rf(c 5 ) should be relat- 
ed as 

T](c s )=(?p ul (c s ) . (21) 

This result has been experimentally verified with inor- 
ganic and organic aqueous solutions such as NaCl, KC1, 
urea, ADP (NH 4 H 2 P0 4 ), KDP (KH 2 P0 4 ), TGS 
[(C 3 H 5 N0 2 ) 3 H 2 S0 4 ], and glycine [4-11] taken at 25 °C 
and normal pressure. For all these solutions it was found 
that the dependence of their bulk densities pic {X> ) 
on the bulk solute concentration was linear: 

pic ^ )=a 0 + a ] c x (Table I gives coefficients a 0 and a x 
for different solutions). The error of such a linear inter- 
polation of the experimental density data was always 
within 0.01%. The experimentally obtained dependences 
of the solution shear viscosity r\ic ^ ) on the bulk solute 
concentration c „ for the TGS and glycine aqueous solu- 
tion are presented in Figs. 1(a) and 1(b) (the experimental 
error of the viscosity measurements is within 15%). In 
these figures the solid lines correspond to the viscosity ex- 
perimental data versus solute concentration whereas 
short-dashed lines represent the sample function 
^7 sample^ <x> ) — (?p W2 ic ^ ). It follows from the straightfor- 
ward comparison between the viscosity experimental data 
line and the sample function line that their intersection 
approximately corresponds (the error in correspondence 
is within 15%) to the saturation concentration c s at the 
given temperature and pressure for every tested solution 
(see Table 1). Therefore, the analytically derived con- 
clusion that at the vicinity of the saturation point there is 
the specific relationship, given by expression (21), be- 
tween the solution viscosity and the density is experimen- 
tally confirmed with an accuracy of 85%. 


1.35 


1.25 


saturation point 


■ ftl! ' li i 

20 % 


V. SOLUTE FLOW TOWARDS 
THE GROWING CRYSTAL SURFACE 


1.15 \- 
1.05 ■ 

0.95 ! 

f 

0.85 J ■ 

0 0.05 


solute concentration, c (in mass fraction) 


0.1 0.15 C 0.2 

ut 


FIG. 1. Dependence of the bulk solution viscosity r\ (c w ) (cP) 
on the bulk solute concentration c x (mass fraction) for the TGS 
and glycine aqueous solutions. Solid lines correspond to experi- 
mental data, whereas short-dashed lines correspond to interpo- 
lation by means for the sample function r/ sarnp , c (c w ). 


To define a complete system of equations describing 
isothermal solute diffusion in the natural convection case 
one has to supplement the general Navier-Stokes and 
continuity equations (4a) and (4b) by the corresponding 
solute diffusion equation. For the particular case where a 
supersaturated solution can be considered as an in- 
compressible fluid and in the stationary limit, the two- 
dimensional equation for the convective solute diffusion 
acquires the form 


3c , 3c 3 

* U 1 + ^ V 2 = ^ 

dx [ 3x 2 3x | 


D(c) 


8c 

dx, 



Die) 


8c 
8 a 2 


( 22 ) 
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It is assumed in this equation that for the saturated and 
the supersaturated solutions the solute mass diffusivity 
Die) is dependent on the solute concentration c. The evi- 
dent analytical form of this dependence is not established 
yet. However, numerous experimental investigations 
[4-11 ] have demonstrated that D(c) is a nontrivial and 
strong function of the solute concentration c. In order to 
solve Eq. (22) within the DBL, it is assumed that there 
exists an expansion (3) for the diffusivity D(c). The fol- 
lowing utilization of the fact that the DBL thickness 
6^ x { ) is very small compared to the characteristic 
length L of a crystal plate allows one to considerably sim- 
plify Eq. (22): 


3c A 3c 
— v \ + ^r v 2 


3xj 


3x ? 


= D(c„) 


3 2 c 

3c 


3x 2 


0(b dlff (x { )-x 2 ) 


(23) 


where 


_ 3 ln[P (c)] 
3c 


It is noteworthy that y „ < 0 since in the supersaturated 
region the solute mass diffusivity Die ^ ) is a decreasing 
function of the solute concentration c^: the diffusivity 
Dlc^ ) declines to zero when c ^ is approaching spinodal 
concentration at the given temperature and pressure. 

To solve Eq. (23) let us introduce the new couple 
(x t ,z) of independent variables, replacing (xj,x 2 ) by 
(X\,z): 


3 

3x, 


3x, 


4x 1 3z 


dx-> 


-Sd r 

Cj ** 3z • 


The following replacement of the velocity components v i 
and v 2 by their expressions given by relations (8) and (9) 
allows one to rewrite Eq. (23) in terms of the new couple 
of independent variables: 


f'(z) 


3c 3/(z) 8c 
4x , 3z 


3xj 


D(c„) 




3 2 c , 

3 c 

2 

3z 2 Yx 

3 z 

0(z M -z) 


(24) 


where the function f(z) within and beyond the DBL is 
given by relations (17a) and (17b), respectively. Taking 
into account the conclusion, obtained in the analysis of 
Eq. (10), that within the DBL the solute concentration c 
is a function of the only variable z, c = c(z), provides the 
following simplification of Eq. (24): 


3 fl .3c _ 


e, 

e 


■D(c.) 


(25) 


In terms of the variable z, boundary conditions for this 
equation can be reduced to the form 


c(z)| z=0 =c(0). 


z = 0 


=0 o [c(O)-c 5 ] , 


(26) 


where c(0) is the solute concentration on the growing 
crystal surface (z =0). It is natural to assume that in the 
crystal growing regime there exists the following double 
inequality: c s 5c(0)<c 0O . The constant j3 0 =/3(c(0)) > 0 
is the c(0)-dependent coefficient that characterizes the 
rate of solute exchange between the crystal surface and 
the solution. The second boundary condition (26) de- 
scribes such a situation where the solute mass flow to- 
wards a crystal surface is positive if c ( 0 ) > c s . 

The solution of Eq. (25) subjected to boundary condi- 
tions (26) is straightforward: 


c(z) = c w (z)0(z oo -z) + c b (z)6(z-z (X) ) . 


(27) 


In this expression c w (z) and c b (z) are the solutions for the 
solute concentration profiles within and beyond the DBL, 
respectively, 

1 


c,Az ) = c(0) + - 


-In 


i+r 


. z/z 


X f dx e 

J n 


(9/4 )Scx * 


cAz)~- 


where Sc — v^/Dic ^ ) is the Schmidt number. The 
analysis of the solution c(z) for z=0 and z^ provides a 
possibility to determine the coefficient /? 0 as the following 
monotonically decreasing (without local extremum) func- 
tion of c(0): 


Po- 


1 

r oo^c 


Mao ~C(0)]__ i , t 

c(0 >_ — |J> 


. (9/4)Scx 4 


(28) 


Substituting this expression for the coefficient (3 0 back 
into relation (27) for the solute concentration c w (z ) we 
obtain 


c u ,(z) = c(0H — — In 
Y oc 


\+(e 


X 




- 1 ) 


X 


Z/Z„ 


dx e 


(9/4)Scjc 


f dx 

J n 


,(9/4)Scx 4 


(29) 


Differentiating this expression with respect to z and hav- 
ing in mind the relationship (1 lb), one can derive the fol- 
lowing condition imposed on the problem charcteristics: 

y oo [c(°)-c aci ] 


F(\)=\~e 
where 


(30) 


3 2 c 

3c 

2 




3z 

9(z„ 0 -z) 

F 

z 


y ~ z/z 

— e” (9/4 ,sc f dxe 


(9/4)Scjc* 


It is noteworthy that F(z/z o0 )< 0 since y ^ <0 in the su- 
persaturated region. Utilization of condition (30) allows 
one to rewrite the result (29) for the solute concentration 
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c w (z) profile as the following function of the variable z: 


c w (z) = c(OH — —In 

Y oc 


1 + 



1 —F( 1 ) 


Ola) 


Taking into account that z =z oc x 2 /6 di jf(x 1 ), one can 
represent expression Ola) for c w (z) in terms of the initial 
coordinate variables x { and x 2 , 


c (z)=c(0H — —In 

7 oo 


1 +- 


1 ~F( 1 ) 


01b) 


Therefore, the most favorable regime for the solute mass 
transfer corresponds to such solute concentration regions 
where the ratio of d \n[D (c)]/dc to d \n[rj(c)]/dc ac- 
quires a maximum value. Expression (34b) for the ratio 
M oo, Sc) can be considerably simplified if one makes use 
of the Watson lemma [15] in the calculation of the in- 
tegral involved. It is straightforward to demonstrate that 
in the limiting case where Sc »f there exists an approxi- 
mation 

f * £ (9/4)Sc 

^ o 9 Sc 

Substitution of this result back into relation (34b) pro- 
vides the following simple expression for the ratio 

M 00 , Sc): 


In the crystal growth problem it is essential to know 
the solute diffusional flux y diff (xi,x 2 ) directed towards 
the growing crystal surface. This flux is defined as 


D(C W ) 


dcjz) 
3 z 


= <? 2 ] X 


V*D(c w ) 


dcjz) 

dx 2 


Let us find separately the flux y diff (z) on the DBL edge 
[x 2 =8 diff (x 1 ) or z=z oc ] and on the crystal surface 
(x 2 =0 or z =0): 


;diff< Jc uS diff ( Jc i)) = - D ^oo )r M , (32a) 

D(c( 0)) e ~ {9/4)Sc 

7diff( * l > 0 ) =./diff < x i , S d iff< * l )) ™ \ \-F(\) ’ (32b> 


where Z)(c(0)) is the solute mass diffusivity on crystal 
surface and (S 2 is given by expression (18). Having in 
mind that within the DBL the solute diffusivity D (c) is a 
weak function of the solute concentration c [an assump- 
tion already used in expansion (3)], one may represent 
£>(c(0)) in the form: 


D(c(0)) = D(O 



Therefore, the ratio Wy^K^Scl^UpO)/ 
y diff [x )], which characterizes the efficiency of the 

solute mass transfer towards the growing crystal surface, 
is given by the expression 


^(r 0 o»^oo»Sc) = We,Sc) = - 


il+€)e 


— (9/4)Sc 


1+6C 


- ( 9 /4 )Sc 


f'd£e ( 9 / 4,Sc t 4 

J n 


(33) 


where e= \ y ^ \ /k^. The analysis of this expression gives 
that the ratio Me, Sc) is the monotonic function of the 
both variables e and Sc. However, Me, Sc) is the increas- 
ing function of the variable e (0 5e5 oo ), whereas it is 
the decreasing function of the variable Sc (0 5 Sc 5 oo ). 
For the given Schmidt number Sc the ratio Me, Sc) ac- 
quires a minimum value when e = 0 and a maximum 
when e— oo : 


min[X(x,Sc)| Sc=consl ] = X(0,Sc) = e (9/4)Sc , 
max[A.(jt,Sc)| Sc=const ] 

= A( oo,Sc)= 


(34a) 

(34b) 


Moo,Sc) = 9Sce“ (9/4)Sc , (35) 

where Sc »j. 

It is the usual situation when the solute mass diffusivity 
Die) is a stronger function of the solute concentration c 
than the solution shear viscosity 77 (c). Therefore, the 
case of e=0 corresponds to a complete disregard of the c 
dependence in diffusivity as well as in viscosity. On the 
contrary, the case of e= 00 describes such concentration 
regions where any change of solute concentration leads to 
the much faster changing of the c-dependent diffusivity 
than of the c-dependent viscosity. It is noteworthy that 
accounting for the c dependence of the solute mass 
diffusivity D{c) and the solution shear viscosity 77 (c) 
leads to an increase of the ratio Me, Sc) estimations 
MO, Sc) < M 00 , Sc). 


VI. CONCLUSIONS 

It follows from the analysis presented above that (a) in 
the vicinity of the saturation concentration c s there exists 
the relationship between the solution shear viscosity and 
density, rj(p s ) a: p ]/2 {c s ); (b) the most favorable regimes 
for the solute mass transfer towards the growing crystal 
surface can be achieved when e»l; and (c) the estima- 
tions of solute mass flow towards the growing crystal sur- 
face performed within the approaches disregarding the c 
dependence of D(c) and 77 (c) lead to its underestimation. 

It is well known that under microgravity conditions 
one may expect a significant improvement in crystal 
growth since the DBL thickness increases with the de- 
crease of the gravity acceleration constant g [see expres- 
sion (19)]. However, as it follows from expression (19), it 
is not necessarily the case. For example, at low supersa- 
turations the ratio (2k ^ — 1 )/g can still be small 

even at microgravity conditions. This prevents the for- 
mation of the appropriate boundary layer needed for im- 
provement of the crystal growth process. Thus, to 
achieve such an improvement of crystal growth one has 
to obtain supersaturations so that 2k ^ /co» 1 . 
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